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IONOSPHERIC  IRREGULARITY  PHYSICS  MODELLING 


I.  Introduction 


The  term  ionospheric  Irregularities  has  become  synonymous  with  the 
terms  plasma  (or  electron)  density  fluctuations,  structure,  and 

striations.  Under  a  wide  variety  of  different  conditions,  ionospheric 
irregularities  are  exhibited  at  all  geomagnetic  latitudes,  longitudes,  and 
almost  all  altitudes.  These  l rregulart t ies  are  generally  magnetic  field 
aligned,  i.e.,  there  is  hardly  any  variation  along  the  geomagnetic  field. 
Thus,  the  irregularity  wavelength  perpendicular  to  the  geomagnetic  field  is 
much  smaller  than  its  parallel  wavelength.  The  two  dimensional  nature  of 
the  irregularities  means  that  detailed  two  dimensional  (perpendicular  to 
the  ambient  magnetic  field)  numerical  simulations  can  he  performed.  There 
is  a  diversity  of  scale  sizes  associated  with  ionospheric  irregularities, 
spanning  some  5-6  orders  of  magnitude  (10's  of  centimeters  to  manv 

kilometers).  ’'ihile  this  diversity  makes  for  a  richness  of  study,  numerical 
simulations,  with  todays  computer  storage,  etc.,  cannot  possibly  handle 
this  broad  scale  spectrum,  i.e.,  the  problem  must  be  examined  in  pieces, 
"evertheless ,  the  inadequacy  of  linear  instability  theories  to  describe 
most  aspects  of  ionospheric  irregularities  and  their  attendant  effects  has 
led  to  studying  from  a  first  principles  physics  point  of  view,  by 

theoretical  and  numerical  simulation  techniques,  the  nonlinear  evolution  of 
the  various  irregularities.  The  study  of  ionospheric  irregularities 
requires:  (i)  a  description  of  the  zeroth  order  (laminar  or  mean  flow) 

state  of  the  ionosphere;  (ii)  a  relevant  instability  mechanism;  and  (iii) 
the  ability  to  follow  the  nonlinear  development  of  such  instabilities. 

The  practical  importance  of  such  studies  of  ionospheric  irregularities 
derives  primarily  from  the  fact  that  such  irregularities  can  degrade 
satellite  C  I  systems  through  scintillation  effects  (e.g.,  for  the  case  of 
F  region  irregularities,  irregularities  with  wavelengths 

~  IOC's  of  meters  to  ~  kilometer  can  cause  scintillation  effects).  The 
studies  should  he  able  to  provide  a  predictive  capability  in  terms  of  where 
these  irregular  structures  can  be  expected  and  what  their  properties  will 
be . 

This  paper  presents  an  overview  of  studies  carried  out,  using 

theoretical  and  numerical  simulation  techniques,  for  ionospheric  F  region 
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irregularities  caused  by  plasma  instabilities.  Results  will  be  presented 
from  studying:  (A)  ionospheric  F  region  plasma  cloud  striation  phenomena 
(generated  by  the  _E  x  gradient  drift  instability);  (B)  equatorial  spread 
F  phenomena  (generated  by  the  collisional  Rayleigh-Taylor  instability);  and 
(C)  high  latitude  diffuse  auroral  F  region  irregularity  phenomena 
(generated  by  the  current  convective  instability  and  F_  x  R_  gradient  drift 
instability).  For  further  study  the  interested  reader  is  referred  to 
references  1-4  for  recent  reviews  of  the  field. 

II.  Theory 

Initially,  we  will  be  interested  in  long  wavelengths  (i.e.,  much 
greater  than  the  ion  gyroradius).  Thus,  the  starting  point  is  the  basic 
plasma  fluid  equations,  which  will  be  applicable  to  all  the  ionospheric  F 
region  irregularity  phenomena  which  we  present.  When  we  study  the 
individual  irregularity  phenomena  stated  in  section  I,  we  will  make  the 
appropriate  approximations.  The  two  fluid  equations  which  describe  the 
electron  and  ion  dynamics  are  as  follows. 
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In  the  above  equations  the  subscript  a  denotes  species  (subscript  e  is 
electron,  l  is  ion),  n  is  density,  _V  is  fluid  velocity,  P  is  production, 
v  is  the  chemical  recombination  coefficient,  p  is  the  scalar  pressure  (p  = 
nT,  T  expressed  in  energy  units),  V  is  the  gradient  operator,  B  is  the 
ambient  magnetic  field  (taken  to  be  uniform),  e^  is  the  charge,  c  is  the 
speed  of  light,  m  is  mass,  _g^  is  gravity,  van  is  the  charged  particle- 
neutral  collision  frequency,  F  is  electric  field,  U  is  neutral  wind 
velocitv,  v  .  is  the  Coulomb  collision  frequency,  and  J  is  current. 
Kquation  (4)  can  be  obtained  from  (1)  by  using  the  appropriate  form  of  (1) 
for  electrons  and  ions,  subtracting  the  two  equations  and  setting 
ne  a  nl  °  n*  fqnat ion  (1)  is  the  continuity  equation,  (2)  is  the  momentum 
equation,  and  (4)  is  the  divergence  of  the  current. 

We  briefly  outline  the  mathematical  representations  (solution  methods) 
that  are  involved  and  basically  common  to  all  three  processes  (plasma  cloud 
striations,  equatorial  spread  F,  and  high  latitude  diffuse  auroral  F  region 
Irregularities).  We  use  numerical  simulations  of  the  appropriate  system  of 
coupled  nonlinear  differential  equations  which  are  derived  theoretically 
and  thought  to  represent  the  process  being  studied  (variations  of  (1)- 
(5)).  Starting  with  the  three  dimensional  continuity  and  momentum 
equations  for  the  ions  and  electrons  (1),  we  note  that  the  spatial  and 
temporal  scales  of  physical.  Interest  are  as  such  to  allow  us  to  make  the 
following  simplifications.  The  electric  field  is  elect ros tat ic  and  thus 
representable  hy  F=  — V 4» .  The  geometry  is  basically  two-dimensional,  t.e., 
the  negligible  components  of  7$  and  the  ion  velocity  parallel  to  _B  allow  us 
to  reduce  the  ID  problem  to  a  series  of  20  problems  coupled  hy  the 

invariance  (or  near  invariance)  of  <(>  along  a  field  line  (to  he  sure  for  the 

high  latitude  diffuse  auroral  irregularities  with  the  current  convective 
instability  special  care  must  he  taken).  Ouasi-neutrality  holds,  t.e., 

ne  81  nt  ,ln^  results  in  (4).  4fter  making  the  above  simplifications  and 
neglecting  inertial  effects  we  are  left  with  two  general  types  of  two 
dimensional  partial  differential  equations:  (a)  one  or  more  nonlinear 

continuity  equations  for  ions  or  electrons  (hyperbolic,  initial  boundary 

value)  and  (h)  a  single  linear  equation  for  the  electrostatic  potential 

(elliptic,  boundary  value,  large  sparse  matrix)  which  is  a  Poisson-like 
equation  resulting  from  V*J  =  0. 
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III.  Examples  of  Ionospheric 
Irregularities  Studied 


A.  Plasma  Cloud  Strlations 

Experimental  studies '’“^of  artificially  Injected  plasma  clouds  Into  the 
ionosphere  have  provided  a  great  amount  of  information  concerning  not  only 
ambient  ionospheric  conditions,  e.g.,  electric  and  magnetic  fields,  but 
also  the  structure  and  morphology  of  evolving  plasma  clouds  themselves  by 
means  of,  for  example,  scintillation  and  power  spectrum  studies.  The 
characteristic  plasma  cloud  initial  steepening,  elongation,  and  striatlon 
formation  have  been  explained  by  applying  the  linear  theory  of  the  E  x  _B_ 

O 

gradient  drift  instability,  originally  developed  for  laboratory  plasmas  , 
to  plasma  cloud  geometries.  More  recently  the  nonlinear  stabilization 

of  the  long  wavelength  E_  x  gradient  drift  instability  in  ionospheric 
plasma  clouds  has  been  studied.^  Numerical  simulation  studies^-^!  of 
barium  plasma  clouds  with  a  background  ionosphere  have  reproduced  not  only 
many  of  the  gross  observational  features  of  plasma  cloud  evolution,  but 

99  9"\— oa 

also  their  spatial  power  spectra  ,  minimum  scale  size4  4  and  outer  scale 

9  S 

size  or  correlation  length.  J  Furthermore,  numerical  simulations  of  the 
local  E_  x  B_  gradient  drift  instability  in  ionospheric  plasma  clouds4  have 
yielded  spatial  power  spectra  and  saturated  wave  amplitudes  that  are 
consistent  with  experiment. 

Figure  1  depicts  the  basic  _F_  x  _B_  gradient  drift  instability  geometry 
for  studying  plasma  cloud  strlations.  We  are  showing  a  simplified  one 
dimensional  geometry  for  the  initial  plasma  cloud  density  profile  (or 
Pedersen  conductivity),  N(y).  In  reality  the  cloud  is  initially  two 
dimensional  perpendicular  to  the  ambient  magnetic  field,  _B_.  However,  for 
simplicity  and  the  fact  that  the  initial  one  dimensional  cloud  geometry 
develops  two  dimensional  (1  _R)  structure  (strlations)  similar  to  the 
initial  20  cloud,  we  have  chosen  to  use,  in  figures  1  and  2,  an  initial  11) 
plasma  cloud.  In  figure  1  the  electric  field,  F_,  is  in  the  x  direction  so 
that  x  _B  is  in  the  negative  y  direction  and  the  so-called  back  side  of 
the  cloud  is  on  the  positive  y  side.  This  striatlon  geometry  is  equivalent 
to  having  a  neutral  wind,  U,  blowing  in  the  y  direction  and  having  no  _F. 
In  figure  1  k  represents  a  perturbation  vector. 
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j  (C)  (d) 

j  Fig.  2  —  Heal  space  isodensity  contour  plots  of  !(x,y)/I()  for  L  =  (>  km  using  the  d'J  random  ini¬ 

tial  conditions  at  (a)  t  =  0  sec.  (I>)  t  =  200  sec,  (c)  t  =  1000  sec.  and  (d)  t  2400  sec  (after  ref.  20). 


For  barium  clouds  released  at  200  km  or  more  v  /Q  <<  1, 

a  a 

where  H  (eaB/mac).  The  electrical  conductivity  of  the  plasma  is 

dominated  by  the  Pedersen  conductivity, *  (nec/TOCv^/Hj).  Due  to  the 

very  high  conductivity,  a.  ,  along  magnetic  field  lines 

6  ^ 

(typically  ~  10  )  ,  the  magnetic  field  lines  may  be  regarded  as 

equipotentials  and,  as  a  result,  a  field  line  integrated  (2D)  model  is 
justified.1*-^  Furthermore,  for  large  plasma  clouds  (i.e.,  large  magnetic 
field  line  integrated  Pedersen  conductivity  compared  with  that  of  the 
background  ionosphere)  the  cloud  interaction  with  the  background  level 
(second  level,  i.e.,  a  separate  continuity  equation  and  momentum  equation 
for  ionospheric  ions)  is  ignored ^ ^ * »23-24.  Thus,  using  an  x,y,z 
coordinate  system  with  _B  =  B^  (B  also  uniform),  ambient  electric  field 
and  after  transforming  to  a  frame  drifting  with  the  x  B^  velocity,  = 
(cRq/B)  the  2D  model  equations  for  the  magnetic  field  line  integrated 
plasma  cloud  Pedersen  conductivity  E  (x,y)  (where  E(x,y)  =  fa  dz)  and  the 
self-consistent  plasma  cloud  electrostatic  potential  <j>,  the  set  of 
equations  (l)-(5))  can  be  written  as 1 2 , 2 1 , 23-24 

|f-  +  |  (£  x  V$)*7E  =  DV2E  (6) 


where  V<j>  =  --E_( x , y > ,  Fj_(x,y)  is  the  total  electrostatic  field, 

V  =  (3/3x,  3/3y),  and  D  is  the  cross  field  diffusion  coefficient.  For 
barium  plasma  clouds  D  lies  in  the  range  0.6  -  6  m  /s.  Fquation  (6) 
results  from  the  magnetic  field-line  integration  of  the  ion  continuity 
equation  while  (7)  is  derived  from  the  current  conservation,  (4).  Also  in 
obtaining  (6)  and  (7)  the  neutral  wind  has  been  set  to  zero.  If  we  set  0  = 
0  then  we  have  the  simplest  set  of  equations  describing  striation  phenomena 

E  =  E  +  E, 
o 


in  F  region  ionospheric  plasma  clouds. 

By  linearizing  (6)  and  (7)  i.e., 


etc.,  and 


assuming  E,  $  *  exp[i(k  y  +  k  x  -  wt)]  with  k*R  =  0,  kl.  »  1  we  find 

y  x  - 

usual  F  x  B  gradient  drift  instability  growth  rate 


cF  k  2 

v  = _ 2.  f_y.i  _  Dk^ 

Y  BL  k 


where  k2  =  k2  +  k2,u)  =  o)  +  iy,  L-^  =  3  In  E  /3x.  Fquation  (9)  shows 
x  y  r  *  o 

that  only  one  side  of  the  plasma  cloud  is  unstable  (positive  L) . 

Numerical  Simulations.  F.quations  (6)  and  (7)  were  solved  numerically 
on  a  mesh  consisting  of  258  grid  points  in  the  x  direction  (J^  x  _B_ 
direction)  and  102  points  in  the  y  direction.  With  a  constant  grid  spacing 
of  310  meters,  the  real  space  dimensions  of  the  mesh  were  BO  km  along  x  and 
31  km  along  y.  The  magnetic  field  line  integrated  Pedersen 

conductivity  E(x,y)  in  (6)  was  advanced  in  time  using  a  multi-dimensional 
flux-corrected  variable  time  step  leapfrog-trapezoid  scheme  which  is 
second  order  in  time  and  fourth  order  in  space.  At  each  timestep  the  self- 
consistent  electrostatic  potential  $  due  to  the  ion  cloud  was  determined 

90  9Q 

from  (7)  using  a  Ohehychev  Iterative  method  •  with  a  convergence 
criterion  of  10-^.  Periodic  boundary  conditions  were  imposed  along  the  y 
direction  with  Neumann  conditions  (3/3x  =  0)  along  the  x  direction.  These 


boundary  conditions  result  in  a  realistic  representation  of  plasma  inflow- 
outflow  at  the  boundaries  in  the  x  R_ direction. 

Initially,  the  field  line  integrated  Pedersen  conductivity  of  the 
plasma  cloud  was  taken  to  be  of  the  form 


E(x,y,t=0)/E  =[M  exp(-x2/Lz)+D. I] ( l+t(x,y)) 
o 


Numerical  studies  have  been  done  where  e(x,v)  was  given  an  rms  value  of  V 

93-25 

and  generated  from  a  randomly  phased  Gaussian  power  spectrum  .  Also, 

<1  A 

two  other  cases  have  been  studied  :  fn  one  case  e(x,y)  =  A  cos  3kpV,  a 

single  monochromatic  wave  along  the  linearly  most  unstable  y  direction  with 

wavelength  A  =  2ir/3k  =  in  km  where  k  =  2rr/3n  km  *  and  the  other 

F  F 

case  e(x,v)  =  A(l-2r(x,v))  where  r(x,v)  is  a  random  number  between  0  and 
1.  The  second  case  models  the  many  wave  initial  condition  (it  produces 
an  e(k)  which  Is  essentially  flat  in  k  space)  and  it  is  this  one  which  we 
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present  here. 

Figures  2a-2d  illustrate  the  evolution  of  an  initially  slab-like 
barium  plasma  cloud  driven  unstable  hv  using  purely  random  initial 
conditions  with  maximum  amplitude  of  17  (A  =  0.03)  for  L  =  6  km  (Reference 
30  also  presents  an  A  =  0 . 1 S 3 .  Also  for  this  simulation  VQ  =  100  m/s,  0  = 
Im'Vs,  and  M  (in  eqn.  (10))  =  1,  so  that  the  maximum  ratio  of  the  field 
line  integrated  cloud  Pedersen  conductivltv  to  the  background  ionosphere  is 
approximately  10.  In  this  simulation  the  initial  evolution  of  a  test  wave 
is  influenced  not  only  by  the  ambient  plasma  cloud,  but  also  by  a  many  wave 
background.  Figure  2b  displays  the  cloud  at  t  =  200  sec  and  shows  the 

initial  random  perturbations  developing  on  the  backside.  Figure  2c  shows 
the  field  line  integrated  Pedersen  conductivity  contours  at  t  =  1000  sec 
with  striation  development  and  elongation  in  evidence.  Figure  2d  shows  the 
cloud  configuration  at  t  =  2400  sec  and  is  similar  to  the  late  time 

configurations  using  other’  initial  perturbation  conditions. 

Figure  3  shows  sample  ID  power  spectra  of  plasma  cloud  Pedersen 

conductivity  both  parallel  P(l<x)  and  perpendicular  P(k^)  to  the  x  _B_ 

drift  corresponding  to  t  =  2400  sec  of  figure  2.  These  power  spectra  were 

obtained  by  first  Fourier  transforming  the  real  space  cloud  Pedersen 

conduct  ivi  tv  6£(x,v)  -*•  SE(k  ,  k  )  where  <SE(x,v)  =  E(x,v)-E  ,  T  the  maximum 

x  v  0  0 

cloud  conductivity.  The  power  spectral  density  1 6E(k^,kv)/E^|2  was  then 
formed  and  the  one-dimensional  power  spectra  P(kx)and  P(kv)  were  computed 
whe  re 


P(k  )  =  /dk  |4£(k  ,k  )/£  |2 
x  y  x  y  o 


(11) 


and 


P(k  )  =  /dk  | 6£(k  ,k  )/E  |2  (12) 

v  x  x  v  o 

These 
with 
scale 


transverse  averaged  power  spectra  P(kx)  and  P(ky)  were  then  fitted" 

a  three  parameter  (spectral  strength  P  ,  spectral  index  n  and  outer 

wavenumber  k  )  power  law  of  the  form 
oor 
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(13) 


p(k  )  =  P  (1  +  k  /k  ) 
a  oa  a  oa 


- n  /  2 
a 


where  a  =  x  or  y.  The  method  used  to  extract  the  best  fit 
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parameters  V  ,  n  ,  k  is  a  nonlinear  least  squares  procedure  which 

oa  a’  oa  M  r 


computes  P  and  n  directly  and  then  iterates  to  find  k 

oa  a  J  oa 


The  time 


histories  of  the  best  fit  spectral  indices  n  and  nv  are  displayed  in 

x  y 

figure  4.  After  initial  transients,  the  spectral  index  in  the  F^  x  B 
direction  becomes  n  3  2  with  n,  3  2-3.  These  spectral  indices  are  in 

x  y 

7  ^  1 

agreement  with  both  experimental  values  ’  and  previous  one 

oa  pc  29 

level  *  and  two  levelz/'  numerical  simulations  of  ionospheric  barium 

clouds  using  different  initial  conditions. 

Summary.  Plasma  fluid  equations  have  been  used  to  studv,  by  numerical 

simulation  methods,  .the  nonlinear  evolution  of  plasma  cloud  striation 

phenomena,  where  the  striations  are  driven  by  the  E  x  B  gradient  drift 

instability.  The  one  dimensional  power  spectra  shown  in  figures  3  and  4 

exhibit  a  power  law  spectrum,  i.e.,  the  x  power  spectrum  is  « 

with  n  “2  for  2n/k  between  1  and  80  km,  while  the  y  spectrum  «  k  ny 
xx  ’  J  r  y 

with  ny  3  2-3  for  2m/k  between  1  and  10  km.  The  numerical  results  are 
consistent  with  recent  experimental  *  ,  theoretical  »  and  other 

numerical  simulation  '  studies  of  ionospheric  plasma  clouds  (reference 
26  discusses  striation  wave-lengths  in  the  100m  -  1  km  regime).  Striation 
studies  to  resolve  shorter  and  shorter  wave-lengths,  as  well  as  the 
stopping  ^f  striation  bifurcation  and  coupling  to  other  ionospheric  levels 
are  ongoing. 


B.  Equatorial  Spread  F 

Equatorial  spread  E  (ESF)  was  discovered  over  four  decades  ago^  as 
diffuse  echoes  on  ionograms.  Notwithstanding  this  fact,  it  is  only 
recently  that  our  understanding  of  this  phenomena  has  increased 
dramatically  (see  references  1-4).  Before  we  proceed  let  us  examine  the 
basic  equatorial  F  region  geometry.  Figure  5  depicts  the  geometry  under 
which  ESF  occurs.  N(y)  represents  the  background  zero  order  electron 
density  as  a  function  of  altitude  (y).  Gravity  _g_  points  down,  the  ambient 
geomagnetic  field  B_  is  horizontal  and  pointing  north  (z),  and  k  represents 
a  horizontal  perturbation  vector  and  points  westward  (x).  The  maximum  in 
electron  density  defines  the  F  peak.  The  bottomside  of  the  profile 
steepens  at  night  due  to  chemical  recombination  effects  and  upward  ambient 
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F  x  B  electrodynamic  forces.  One  might  expect  this  equatorial  Ionospheric 
geometry  to  he  unstable  to  a  variety  of  plasma  Instabilities  because  in  the 
plasma  physics  community  Figure  5  represents  the  classical  flute  mode 
geometry.  It  Is  reminiscent  of  a  heavy  fluid  (higher  density  region  on  the 
hottomside)  being  supported  by  a  light  fluid  (low  density  region  on  the 
bottomside). 

In  the  past  few  years  much  effort  has  been  expended  to  describe  and 
explain  equatorial  spread  F  phenomena  both  from  an  experimental  and 
theoretical  viewpoint*-^.  It  is  now  generally  believed  that  the 
rollisional  Ravleigh-Taylor  instability  mechanism  (or  possiblv  the  F  x  _B_ 
gradient  drift  instability;  in  any  case  a  fluid  type  gradient  instability) 
initiates  equatorial  spread  F  in  the  bottomside  evening  equatorial  F  region 
ionosphere Several  experimental  results  lend  credence  to  this 
idea ^ .  The  generation  of  this  instability  on  the  bottomside  leads  to 
the  formation  of  plasma  density  depletions  (bubbles),  which  then  rise  to 
the  topside.  These  buhbles  were  predicted  by  the  nonlinear  numerical 
simulation  studies  (in  the  collisional  Rayleigh-Taylor  regime,  denoted  CR- 

O  Q 

T)  of  Scarmapieco  and  Ossakow  and  have  been  observed  experimental¬ 
ly.  1  The  initial  nonlinear  numerical  studies^  showed  how  a 

plasma  mode  (CR-T),  which  is  linearly  unstable  on  the  bottomside  only, 
could,  by  nonlinear  polarization  _F  x  _R  forces  in  the  equatorial  geometry, 
produce  bubble  induced  irregularities  on  the  topside,  where  the  mode  is 
linearly  stable.  These  initial  simulation  studies  represented  a  major 
break-through  in  the  undertaking  of  FSF  phenomena  for  they  answered  the 
following  age  old  questions.  (1)  Why  are  plasma  density  fluctuations 
observed  on  hoth  sides  oc  the  F  peak  when  linear  theory  of  the 
representative  instabilities  predict  only  bottomside  growth  (unstable)? 
(?)  Why  does  the  onset  of  topside  F.SF  lag  the  onset  of  hottomside  FSF? 
Other  numerical  simulation  studies  of  large  scale  bubble  phenomena  have 
been  performed^-^4  as  well  as  intermediate  wavelength  studies^  ^  which 
examine  a  chunk  of  the  hottomside  FSF  ionosphere  to  obtain  power  spectra 
for  irregularity  wavelengths  in  the  ln0  m-1  km  regime.  These  studies  help 
to  answer  the  following  additional  questions.  (?)  Mow  are  the  large  plasma 
depletions  (hubbies)  produced?  (4)  '-That  causes  the  westward  tilts  of  the 
bubbles  and  the  observed  radar  haokscatter  plumes?  ( "> )  FTiat  is  t lie  power 
spectrum  for  the  plasma  irregularities  during  FSF? 


1  l 


We  will  now  apply  (l)-(5)  to  the  two  dimensions  perpendicular  to  B^  at 
the  geomagnetic  equator,  making  various  approximations.  If  inertial  terms 
are  neglected  (i.e.,  the  left  hand  side  of  (2)),  take  v  ^  ^  (valid  in 
the  F  region  approximation),  neglect  Ra  in  (2),  neglect  JJ,  and  use  the 
quasi-neutral  approximation,  ne  “  n^  «  n,  then  ( 1 ) — ( 2 )  become  for  the  two 
dimensions  (x,y)  perpendicular  to  _B 
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where  we  have  neglected  neutral  wind  effects,  =  eB/m^c,  _B^  =  B2_  and  we 

have  set  Tg  =  =  0  for  simplicity  (see  reference  42).  It  should  be  noted 

that  the  term  v  n  on  the  right  hand  side  of  (14)  represents  an  effective 
k  ao 

ionization  source  such  that  n  (t)  =  n„  (0),  i.e.,  that  the  zero  order 

ao  ao 

distribution  changes  on  a  much  longer  time  scale  than  any 

perturbations.42’42  The  g_  x  ^motion  in  (15)  has  been  neglected  compared 

with  the  same  ion  motion  because  they  are  in  the  ratio  me/m.£.  Making  the 

electrostatic  approximation  _E  =  —  V d>  and  setting  4>  =  4>0  +  where  the 

subscript  zero  refers  to  equilibrium  or  zero  order  quantities  and  the  tilde 
again  refers  to  perturbed  quantities,  we  obtain  (with  ne  “  K  n) 
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Tn  arriving  at  (18)  we  have  put  (IS)  and  (16)  into  (4)  and  (5)  and 

obtained  V =  m.g/e.  Equation  (9)  can  he  viewed  as  the  electron 

o  l — 

continuity  equation  in  the  -  x  _R  drift  frame  (i.e.,  _g_  x 

9J  frame).  Equation  (18)  results  from  (4).  Equations  (17)  and  (18)  are 

to  he  compared  with  (6)  and  (7).  Tn  their  simplest  forms  (setting  v  =  D  = 

K 

0),  the  equations  are  quite  similar.  Tn  (18)  if  we  think  of  v  as  a 

In 

constant  then  (18)  is  identical  with  (7)  if  we  interpret  (B/c  )  g  x  2  as 

an  effective  electric  field.  To  obtain  a  linear  growth  rate  from  (17)  and 

(18)  we  set  n  =  nQ(y)  +  n  and  put  n,$  a  exp  i(k*x-  o)t ) ,  where  k' ■  =  k  1  + 

9 

ky  andw  r  +  iy.  The  growth  rate  y  then  is,  for  kL  >>  1, 
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where 
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-1 


3  In  n  /3y.  For  horizontal  waves,  i.e.,  k  =  k  this  reduces  to 
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Both  (19)  and  (2D)  show  that  linear  growth  is  possible  only  on  the 
bottomside  of  the  equatorial  ionosphere  depicted  in  fig.  5. 

Two-Dimensional  Olohal  Numerical  Simulations.  Tn  this  section  we  will 
outline  some  two  dimensional  (1  _R_  and  at  the  geomagnetic  equator)  numerical 
simulation  results  of  the  collisional  Rayligh-Taylor  mechanism^.  The 

simulation  plane  encompasses  an  altitude  extent  of  over  200  km  (see  fig. 
6),  from  the  bottomside  to  the  topside  of  the  ESF  region.  Ionospheric 
profiles  of  n  (y),  v,  (v),  and  v  (v)  are  used  in  the  simulations.  A 

in  K 

simulation  mesh  of  140  (altitude,  v)  by  40  (east-west,  x)  was  employed 

with  Ax  =  3  km  (large  case,  f,).  Tn  reference  43  studies  were  also 

performed  with  Ax  =  200  m  (small  case,  S).  The  boundary  conditions  were 

periodic  in  x  and  y  transmi tt ive  for  n  (3n/3v  =  0)  and  Neumann 

(3<fc/3y  =  0)  for  the  induced  potential^.  Tn  these  particular  simulations 

n  was  set  to  zero  in  (17)  so  that  n„  obeved  3n  /3t  =  -  v>  n  and  ''IS)  was 

°  o  -  o  R  o 
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retained  in  tact.  Two  different  types  of  initial  perturbations  and 

different  wavelengths  corresponding  to  the  large  case  and  small  case  were 

used  (four  cases  in  all).  However,  here  we  will  only  present  case  1L, 

which  is  defined  with  Ax  =  5  km  and  a  perturbation  of  the  form  (see  ref.  43 

for  more  details). 


n(x,y,0) 

nQ(y»°) 


0<  |x  |  <  8Ax 


n(x,y,0) 

nQ(y,0) 


1  -e 
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cos 


fiAx  <|xj<16Ax 


n(x,y,0)  . 

nQ(y,o) 


|x | >16Ax 


(21) 


In  contrast  to  this  more  localized  type  of  perturbation  which  has  been  used 
before^  the  pure  sinusoidal  one,  also  used  in  ref.  43,  will  not  be  shown 
here.  [Note:  The  perturbation  for  the  pure  sinusoid  is  given  as 

n(x,y,0)/no(y ,0)  =  1  -  e-^  cos  (ttx/20Ax).  This  different  type  of  non- 

localized  perturbation  was  chosen  to  contrast  with  (21).  The  results 
produced  for  bubble  rise,  etc.  are  similar  in  nature  with  some  slight 
differences  over  the  results  obtained  using  (21).  For  the  sake  of  brevity 
the  reader  is  referred  to  reference  43  for  more  details.) 

Figure  6  shows  isodensity  contours  of  calculation  1L  at  times  300, 
700,  1000  and  1200  sec  after  initialization.  While  the  initial 

perturbation  (21)  is  put  over  the  entire  altitude  (y)  extent  of  the  mesh, 
we  see  from  300  to  700  sec  that  the  perturbation  grows  on  the  bottomside, 
but  damps  on  the  topside  as  it  should  according  to  (20).  The  build-up  of 
the  bubble  is  evident  in  the  1000  and  1200  sec  plots.  At  1200  sec  figure  6 
shows  a  single  plume  of  depleted  ionization  with  no  secondary  central 
bubble  or  side  bubble.  In  figure  7  a  plot  of  n(x,y)/nQ(y)  for  case  lb  at 
1200  sec  is  exhibited.  The  level  of  depletion  is  greater  than  99.9 7-  for 
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Fig.  6  —  Sequence  of  four  plots  showing  isoelectron  density  contours  of  calculation  1L  at  300, 
700,  1000,  and  1200  sec.  Superimposed  on  each  plot  is  a  dashed  line  showing  nQ(y,t)  and  labeled 
at  the  top.  Electron  densities  are  given  in  units  of  cm'3.  The  observer  is  looking  southward  so  that 
B  is  out  of  the  figure  (after  ref.  43). 
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the  entire  10  by  70  km  oval  hole  located  inside  the  tenth  solid  contour  and 
represents  at  least  a  three  orders  of  magnitude  decrease  (bite-out)  in 
plasma  density.  The  bubble  rise  velocity,  computed  from  the  last  two 

frames  of  fig.  b,  is  230  m/s.  These  large  bite-outs  and  large  rise 
velocities  are  in  agreement  with  the  observations  of  McClure  et  al."^  The 
reason  that  the  bubble  is  so  depleted  in  the  1L  case  can  be  seen  with  the 

aid  of  fig.  8.  This  shows  the  induced  electrostatic  potential  at  700 

sec.  We  notice  that  the  fringe  field  surrounding  the  bubble  extends  to  low 
altitudes  (lower  than  say  the  18  case  in  ref.  43).  Since  the  contours  of 
constant  potential  are  in  fact  flow  streamlines,  much  plasma  is  being  drawn 
up  from  low  altitudes  (see  fig.  6).  the  effect  is  similar  to  the  fringe 
field  of  a  parallel  plate  capacitor.  In  the  case  depicted  in  fig.  8  the 
role  of  plate  separation  is  played  by  the  horizontal  scale  size  of  the 

bubble  (see  ref.  43  for  more  details).  We  also  note  that  the  induced 
potential  is  such  as  to  make  the  induced  electric  field  point  from  right 
(west)  to  left  (east)  so  that  the  induced  F,  x  R  motion  is  upward. 

Sunmary .  In  this  section  we  have  shown  how  plasma  fluid  equations  are 
used  to  study,  by  numerical  simulation  methods,  the  nonlinear  evolution  of 
equatorial  spread  F  phenomena,  driven  by  the  collisional  Rayleigh-Taylor 
instability  mechanism.  To  keep  things  simple  we  have  exhibited  the  pure 
collisional  R-T  regime,  i.e.,  horizontal  zero  order  neutral  winds  and/or 
electric  fields  have  not  been  included  in  the  study.  The  results  gleaned 
from  these  numerical  simulations  of  large  horizontal  wavelength  initial 
perturbations  in  the  collisional  R-T  regime  are  as  follows.  The  large 
scale  length  initial  perturbations  evolve  nonlinearly  into  large  horizontal 
scale  length  ESF  bubbles.  These  bubbles  rise  up  from  the  bottomside  to  the 
topside  and  evolve  approximately  on  the  same  time  scale  as  do  the  smaller 
horizontal  scale  length  counterparts^’^.  The  plasma  comprising  these 
large  hubbies  has  its  origin  at  very  low  altitudes  (much  lower  than  that  of 
the  smaller  horizontal  scale  length  bubbles)  and  this  results  in  plasma 
density  depletions  in  the  large  bubbles  very  close  to  1007,  i.e.,  several 
orders  of  magnitude  depletions.  These  results  are  consistent  with 
experimental  results ' .  Studies  to  investigate  shorter 
wavelength  FSF  phenomena,  via  analytic  and  numerical  simulation  methods, 
have  been  ongoing  ,  including  the  decay^  of  these  F.SF 

irregularities.  Large  scale  global  numerical  simulations^  have  been 
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Fig.  8  —  Contours  of  constant  induced  electrostatic  potential  at  700  sec  for  calculation  1L.  Small 
dashed  lines  are  for  negative  potentials  and  solid  lines  are  for  positive  potentials.  The  potential 
drop,  across  adjacent  contours  is  a  constant  except  that  the  zero  potential  is  not  plotted.  The  con¬ 
tours  are  chosen  to  space  the  twelve  contours  (plus  the  zero  contour)  uniformly  from  maximum  to 
minimum  potential.  The  targe  dashed  line  depicts  n0(y)  at  700  sec  (after  ref.  43). 
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started  which  incorporate  the  effects  of  background  Pedersen  conductivity 
(another  level,  similar  to  plasma  cloud  striation  studies)  and  an  F  region 
neutral  wind,  at  the  equator,  blowing  to  the  east.  Results  from  this 
study^  show  that  the  first  effect  can  slow  down  FSF  and  the  attendant 
bubble  rise  phenomena  and  the  neutral  wind  effect  in  combination  with  the 
first  effect  can  result  in  a  westward  tilt  of  the  bubbles  and 

attendant  radar  backscatter  plumes^*”^.  Also,  the  studies^  show  that 
solution  of  the  zero  order  or  equilibrium  equations  result  in  a  plasma 
velocity  shear  as  a  function  of  altitude.  This  shear  is  in  agreement  with 
recent  radar  backscatter  studies^-^  of  FSF  phenomena. 

Fven  though  work  is  still  continuing  on  FSF  phenomena,  a  basic 
physical  picture  of  FSF  that  emerges  through  theoretical,  numerical 
simulation  and  experimental  studies  is  as  follows.  After  sunset  the  F- 
region  begins  to  recombine  and  there  is  effectively  no  F-region 
conductivity  to  short  out  any  irregularities  in  the  F-region.  T\ie  to 
recombination  and  possible  electrodynamic  effects  (upward  motion  of  the  r- 
region  plasma  due  to  _E  x  _B)  the  bottomside  F-region  background  electron 
density  gradient  begins  to  steepen.  The  steepening  due  to  plasma  F_  x  _B 
uplift  caused  by  an  eastward  electric  field  is  much  like  barium  cloud 
steepening  in  equatorial  geometry  (note  a  downward  motion  of  the  neutral 
atmosphere  will  have  the  same  effect).  Uhen  the  altitude  of  the  F-region 
is  high  enough  and/or  the  bottomside  background  electron  densitv  gradient 
steep  enough  to  overcome  recombination  effects,  plasma  density  fluctuations 
will  begin  to  grow  on  the  bottomside  via  the  collisional  Ravi  ei gh-Taylo r 
instability  mechanism  (or  possibly  initiated  by  the  F.  x  B  gradient  drift 
instability).  These  irregularities  will  in  turn  form  plasma  density 
depletions  (bubbles)  on  the  bottomside  which  will  then  nonlinearlv  rise  by 
polarization  F_  x  B  motion  through  the  F  peak  and  cause  topside  spread  f. 
(The  higher  the  F  peak  and/or  the  sharper  the  bottomside  gradient,  and  the 
higher  the  percentage  depletion  the  faster  the  rise  rate  of  the  bubble  and 
the  faster  the  evolution  of  FSF).  These  steepening  bubbles  (on  their 
topside)  can  then  bifurcate  and  form  smaller  scale  structure  hv  a  cascade 
or  two-step  mechanism.  An  eastward  neutral  wind  can  cause  the  huhhles  to 
move  westward  with  respect  to  the  hulk  plasma  motion.  The  long  wavelength 
irregularities  have  a  power  law  PSD  ®  k  ^  '  It  is  present lv  believed  that 
the  small  scale  (<  1  m)  radar  backscatter  observed  irregularities  (kinetic 


type  Instabilities)  are  created  by  the  primary  long  wavelength  fluid  type 
Instabilities  (colltstonal  Ravlpi gh-Taylor,  etc.).  The  turn  on  and  turn 
off  of  these  very  small  scale  Irregularities  and  their  connection  to  the 
longer  wavelength  Irregularities  is  not  well  understood  at  present  and 
requires  further  investigation  (see  ref.  3  for  more  details). 

C.  High  Latitude  Diffuse  Auroral  Irregularities 

DNA  Wideband  satellite  observations  have  exhibited  a  regularly 
occurring  scint 1  Hat  ton  enhancement  in  the  auroral  zone  data,  In  the  region 
of  diffuse  auroral  particle  precipitat  Ion.  The  scintillation  data 

have  indicated  that  the  enhancement  is  due  to  sheet-llke,  L  shell  aligned  F 
region  ionospheric  irregularities^^*1’'.  Tn  the  initial  experimental 
results  ** ' ,  the  irregularities  occurred  where  there  was  a  strong 
northward  total  electron  content  (TFC)  gradient.  Associated  with  the 
northward  TF.C  gradient  is  a  northward  gradient  in  local  F  region  plasma 

density^  as  measured  by  the  Chatanika  radar.  Tn  addition,  measurements 

show  a  d .  c .  electric  field  pointing  westward  or  northwest  and  a  very 

shallow  altitude  plasma  density  gradient, ^  i.e.,  the  northward  TFC 
gradient  dominates.  hater  experimental  results'^  showed  large  scale 
equatorward  convecting  enhancements,  limited  in  latitudinal  extent  with 
steep  poleward  and  equatorward  edges  (two-sided,  much  like  a  large  barium 
c loud  ) . 

At  first  sight,  with  the  dominant  plasma  density  gradient  pointing 

northward,  the  ambient  magnetic  field,  pointing  down  and  the  d.c. 

electric  field  horizontal  (see  fig.  9),  this  would  appear  to  be  a  prime 
geometry  for  the  usual  F  region  F.  x  B  gradient  instability  (see  subsection 
A  and  fig.  1).  however,  with  the  d.c.  electric  field  pointing  westward  or 
northwest,  the  configuration  is  stable  (the  electric  field  would  have  to 
point  eastward  for  instability).  The  saving  feati  •*« ,  however,  is  the 

diffuse  auroral  precipitation  resulting  in  a  current  along  R_,  which  acts  to 
destabilize  the  ahov"  geometrical  configuration  .  The  conditions  of 
having  a  current  along  R_,  a  density  gradient,  Vn  perpendicular  to  can 
result  in  instability  and  these  types  of  instabilities  are  generieallv 
called  current  convective  instabilities'^  (CCT).  The  instability 

has  k(|  <<  (see  next  section)  and  so  the  Irregularities  generated  will  he 
magnetic  field  aligned.  This  current  convective  instability  can  dirertlv 
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result  In  long  wavelength  scintillation  causing  F  region  Ionospheric 
Irregularities  in  the  diffuse  auroral  region.  A  nonlinear  mode  coupling 
theory  has  been  applied  ^  to  the  CCI  to  look  at  saturation  of  the 
instability.  In  addition,  the  linear  theory'’"’  of  the  COT  has  been  extended 

CO 

to  include  magnetic  shear  effects  ,  highly  collisional  (non-magnetized ) 
ions  (e.g.,  the  F  region)'’**,  electromagnetic  effects-’**,  and  ion  inertial 
effects  ,  which  extends  the  instability  to  higher  F  region  altitudes. 
Numerical  simulations  of  diffuse  auroral  F  region  irregularities  have  been 
performed  for  one-sided  density  gradients^®  (where  the  CCI  is  dominant),  as 
well  as  two-sided  density  gradients*^  (where  the  F_  x  B_  gradient  drift 
instability  is  dominant  on  the  poleward  gradient  side,  with  the  CCI  adding 
to  the  instability).  In  this  paper,  in  order  to  emphasize  the  nature  of 
the  CCI  instability,  results  from  the  former  numerical  simulation*’®  will  be 
presented. 

Now  we  present  a  theoretical  model for  a  long  wavelength  fluid  type 
plasma  instability  which  may  account  for  the  scintillation  causing  diffuse 
auroral  F  region  ionosphere  irregularities.  For  our  model  we  take  the 
electron  density  gradient  to  be  pointing  northward  (y),  the  ambient 
electric  field  _E^  is  in  the  westward  direction  (x)  and  the  magnetic  field 
points  downward  (z).  In  our  simple  model  we  have  equated  the  TKC  gradient 
with  a  gradient  in  plasma  density.  We  assume  that  the  horizontal 
(northward)  electron  density  gradient  is  much  sharper  than  the  altitude 
density  gradient,  which  we  neglect.  Our  basic  equations,  from  ( 1 )— ( 5) ,  are 
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and  (5),  where  we  have  neglected  inertial  effects  (left  hand  side  of  (22)), 

2  =  B /  I B | ,  Q  =  |e  |B/m  c,  V  represents  the  diffuse  auroral  particle 

—  —  1  ’  a  a  a  ’  oo  I 
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precipitation  induced  velocity  along  _B_  (  II  and  i  denote  parallel  and 

perpendicular  to  B)  which  results  in  a  zero  order  current  and  the  other 

symbols  have  their  usual  meaning.  In  addition,  we  have  neglected 

temperature  and  neutral  wind  effects  and  included  electron  collisions  to 

first  order  for  completeness.  Equations  (23)  and  (24)  are  valid  for  F 

region  ionospheric  altitudes  such  that  <<  1*  In  (23)  is  taken  to 

mean  ion-neutral  collisions,  while  in  (24)  v  is  really  electron-ion 

e 

collisions.  In  (24)  the  electron  Pedersen  drift  has  been  neglected, 
compared  with  the  ion  Pedersen  drift  in  (23).  Making  the  usual  quasi¬ 
neutrality  assumption,  ne  “  n|  K  n»  final  equations  from  (22),  (23), 

(24)  and  (5)  are 
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(26) 


where  V<J>  =  E  -  with  E(x,y,z,t)  the  total  electric  field  and  = 

V.  „  -  V  „ .  Equations  (25)  and  (26)  are  simplv  a  restatement  of  the  ion- 
ioll  eoll 

continuity  equation  written  in  a  reference  frame  with  velocity  =  - 

E0/B(2  “  v  jO  an<^  quasi-neutralitv  V»  J_  =  n,  respectively.  Tn  (25)  and 
(26)  V  =  (3/3x,  3/3y,  3/3z)  and  =  (3/3x,3/3v).  We  note  the  similarity 
between  (25),  (26)  and  (6),  (7),  and  (17),  (18).  If  we  set  =  0  and 
3/3z  =  0  in  (25)  and  (26)  then  we  recover  the  simplest  form  of  the  F_  x 
gradient  drift  instability,  equations  (simplest  form  of  (6)  and  (7)). 

Equations  (25)  and  (26)  are  then  linearized  such  that  n  =  nQ(y)  +  n,  _E_ 
=  E  8  -V$,  V  =  V  +V  with  the  perturbed  quantities  n, 

V  *  exp  i  [  k  x  +  k  ||  z  -  wt  1 ,  where  u>  =  +  iy.  We  then  obtain 
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where  j  -  veo||*  We  see  that  In  (27)  y  is  independent  of  |k|  and 

only  depends  on  the  angle  that  k  makes  with  _B.  In  the  denominator  of  (27), 

the  first  term  (in  brackets)  multiplying  ^/k^)2  comes  from  the  parallel 

motion  of  the  Ions  and  electrons;  whereas,  the  remainder  of  the  denominator 

comes  from  the  ion  Pedersen  motion.  It  should  be  noted  that  the 

instability  is  essentially  unaffected  by  the  current  direction.  Thus, 

downward  currents  work  just  as  well  as  upward  currents.  In  (27)  if  we  set 

k„  =  0  we  obtain  y  =  -  (n  ^3n  /3y)  (cE  /B)  which  is  the  usual  result  for 
II  o  o  o 

the  E  x  _B  gradient  drift  instability.  For  our  geometry  this  shows  y  is 
negative  which  implies  stability.  For  instability  in  (27)  Y  >  0  which 
implies 
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This  says  that  for  instability,  with  the  westward  ,  we  must  have 

V.k./k  <  D  and  | V , k „ | > ( k  cE  /B)(v. /ft.).  We  can  maximize  this  growth  with 
d  n  x  d  ii  x  it 

respect  to  0  =  k^/k^.  From  (27)  the  growth  rate  maximizes,  in  general  for 
9  =  -  (cF,o/B  Vd)(v1/Sl1)  ± 

{(ceq/b  vd)2(vi/n1)2  + 

(v./n.mn./v  )  +  (0  /v  )1"  }1/2  (21) 
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For  altitudes  ~  350-400  km  corresponding  to  the  observation 

altitudes5^  v^/fi^  ~  10  \  -  10  \  and  E0  ~  10  mV/m  with  ~  -  500 

m/sec  we  find  9  =  9.4  x  10  5.  For  this  case  (27)  gives  a  growth 

-3  -1  -1 

rate  y  -  2.7  x  10  sec  for  L  =  n  (3n  /3y)  ~  50  km.  Including 

max  o  o 

pressure  effects  In  the  problem  Introduces  diffusive  damping  in  (27).  A 
typical  cross-field  diffusion  coefficient,  T)^ ,  Is  ~  0.2  m^/sec  and  a 

Q  O 

parallel  diffusion  coefficient,  D^,  Is  ~  10  m  /sec.  In  the  present  study, 
these  effects  become  important  for  perpendicular  wavelengths,  A^,  <  100  m 
and  parallel  wavelengths,  A  ,  <  1000  km.  However,  typical  scintillation 

"  ~  _4 

causing  perpendicular  wavelengths  are  ~  1  km  and  since  k^/k^  ~  10  we  are 
considering  highly  field  aligned  irregularities.  Larger  parallel  currents, 
due  to  precipitation,  will  of  course  produce  larger  growth  rates.  However, 
too  large  a  current,  Vj  >  1  km/sec,  would  excite  the  collisional  electro- 
static  ion  cyclotron  instability  .  It  may  be  noted  that  the  linear  theory 
of  the  current  convective  instability  proposed  here  favors  a  wavevector 
perpendicular  to  the  TEC  gradient  (as  well  as  B).  However,  in  the 
nonlinear  regime,  mode  coupling  theory  suggests57  that  the  dominant 
wavevectors  will  lie  parallel  to  the  TEC  gradient  (and  local  F  region 
plasma  density  gradient)  thus  accounting  for  the  L-shell  aligned  nature  of 
these  irregularities5^’55.  Basically  in  this  theory,  the  current 
convective  instability  saturates  (stabilizes)  nonlinearly  by  generating 
linearly  damped  harmonics  (i.e.,  those  along  the  northward  plasma  density 
gradient).  These  harmonics  dominate  over  those  between  the  northward  and 

east-west  directions.  In  addition,  these  northward  modes  produce  a  power 

-2 

spectrum  for  the  irregularities  *  k  . 

Diffuse  Auroral  Numerical  Simulations.  In  the  following  simulations 
we  take  advantage  of  the  fact  that  the  fastest  growing,  most  dangerous 
modes  from  linear  theory  are  almost  field-aligned  (k^/k^  <<  1)*  These 
modes  are  of  most  interest  to  us  and,  as  a  result,  we  solve  equations  (25) 
and  (26)  in  a  plane  containing  these  modes  which  is  nearly  perpendicular  to 
the  magnetic  field  while  fixing  the  value  of  kfl/k^  (see  ref.  60).  A 
similar  approach  has  been  adopted  in  numerical  studies  of  drift-wave  and 
trapped-particle^4  instabilities  in  laboratory  plasmas.  The  simulation 
plane  which  is  essentially  horizontal  at  an  altitude  of  350  km  with  a 
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north-south  extent  of  410  km  and  an  east-west  extent  of  160  km  is  identical 


to  the  x'y"*  plane  as  shown  in  Figure  4.  The  system  of  equations  (25)  and 
(26)  was  transformed*^’ ^ *  to  the  x^y'z'  coordinate  system  by  a  simple 
rotation  about  the  y-axis  by  the  angle  6  =  k  ^  /k^  <<  1.  By  solving  the 
equations  (25)  and  (26)  in  the  x'y'z'  system  a  finite  but  small  k^  is 
effectively  introduced  into  the  model.  After  neglecting  the  z"-  dependence 
of  all  quantities,  equations  (25)  and  (26)  were  then  cast  into 
dimensionless  form  and  initialized  with  the  profile  of  the  following 
type  nQ(y")/No  =  (l-A(l-tanh(  y'-yQ)/L))  (1  +  e(x',y'))  where  NQ,  yQ,  and 
L  are  constant  with  A  =  5/11.  This  gives  a  total  density  maximum  to 
minimum  ratio  of  approximately  10.  The  function  e(x'',y')  has  a  root-mean- 
square  value  of  37  and  is  generated  from  a  randomly  phased  Gaussian  power 
spectrum.  The  computational  mesh  consisted  of  258  grid  points  in 
the  y'-  direction  (north-south)  with  102  points  in  the  x'  -  direction  (east- 
west)  so  that  Ay  =  Ax  =  1.6  km.  Periodic  boundary  conditions  were  imposed 

in  the  x“-  direction  with  Neumann  (B/Sy"*  =  0)boundary  conditions  in 
the  y'-  direction.  We  now  drop  the  prime  notation  for  clarity. 

Figure  10  shows  contour  plots  of  n(x,y)/NQ  at  t  =  0,  900,  1400,  1900 

sec.  The  following  set  of  parameters  were  used:  L  =  50  km,  yQ  =  200  km, 

F  =  10  mV/m,  v  /Cl .  =  10  ^  v  /SI  =  10  ^ .  The  value  of  9  =  9.4  x  10  is 

o  i  i  e  e 

held  fixed  and  is  chosen  so  as  to  maximize  the  linear  growth  rate  (27). 
Figure  10  at  t  =  0  shows  the  initial  configuration  which  includes  the  small 
random  perturbation.  Figure  10  at  t  =  900  sec(  y^t  -  5)  illustrates  the 
linear  stages  of  the  simulation  and  shows  unstable  growth  in  the  region 
where  3no/3v  >  0  as  predicted  by  linear  theory.  Figure  10  at  t  =  1400  sec 
exemplifies  the  early  nonlinear  regime  where  lower  density  plasma 

(depletions)  are  moving  in  the  positive  y-direction  (poleward)  while  higher 
density  plasma  (enhancements)  are  convecting  in  the  negative  y-direction 
(equatorward ) .  The  approximate  velocities  of  the  depletions  and 

enhancements  are  270  m/sec  and  30  m/sec,  respectively.  Finally,  well- 
developed  steepened  enhancements  and  depletions  (of  over  907.)  are  seen  in 
Fig.  10  at  t  =  1900  sec.  This  late-time  configuration  is  reminiscent  of 
the  motion  of  depletions  (bubbles)  moving  vertically  in  the  equatorial  F 

90 

region  and  enhancements  (strlations)  in  ionospheric  F  region  plasma 
clouds  .  The  length  scales  in  Figure  10  are  distorted  with  the  depletions 
and  enhancements  longer  and  narrrower  than  is  depicted.  Similar  linear  and 


shuth  ■  Nnn^H 


c  see 


Fig.  10  —  Real  space  isodensity  contour  plots  of  n(x.y)/N0  for  L  -  50  km  at  t  =  0.  900,  1400  and 
1900  sec.  Fight  contours  are  plotted  in  equal  increments  from  1  to  10  with  the  lower  (higher) 
density  contours  denoted  by  solid  (dashed)  lines.  The  magnetic  field  B  is  directed  into  the  figure 
with  the  observer  looking  down  the  magnetic  field  lines  toward  the  earth  (after  ref.  60). 
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nonlinear  development  Is  observed  when  L  =  10  km,  but  on  a  faster  time 

scale.  Figures  1  la-1  lb  give  sample  one-dimensional  spatial  power  spectra 

at  t  =  1000  sec  both  in  the  x-direction  (east-west)  and  in  the  y-direction 

(north-south).  These  spectra  are  obtained  by  first  Fourier 

analyzing  5n(x,v)/N  and  integrating  over  the  direction  in  k-space 

o  — 

corresponding  to  the  north-south  and  east-west  directions,  respectively 
(see  (11)  and  (12)).  For  both  cases,  these  power  spectra  are  well-fitted 
with  an  inverse  power  law.  Similar  power  law  dependences  were  seen  when  L 
=  10  km. 

The  following  physical  picture  of  the  evolution  of  current  convective 

instability  in  the  diffuse  aurora  is  supported  by  these  simulations.  Tn 

the  evening  a  westward  electric  field  _F^  begins  to  form  which  convects 

plasma  in  the  auroral  region  equatorward.  Tn  regions  where  the  northward 

gradient  in  total  electron  content  becomes  well-defined  nearly  field- 

aligned  fluctuations  (kfl/k^  <<  1)  will  grow  unstable  in  regions  where  the 

field-aligned  current  velocities  ,  caused  by  precipitating  particles,  are 

such  that  V,(k11/k1)  >  0  and  |V,k„|>(k.cF  /B)  (v./Sl.).  Tn  the  plane  almost 

dill  d  n  1  o  ii 

perpendicular  to  the  magnetic  field  by  an  angle  9  =  lc (|  /k ^  plasma  depletions 
and  enhancements  will  move  northward  and  southward,  respectively,  while 
steepening  in  the  process. 

Summary.  We  have  investigated  analytically  and  numerically  a  simple 

plasma  fluid  model  to  account  for  the  diffuse  auroral  scintillation  causing 

F  region  ionospheric  irregularities  observed  by  the  T>NA  Wideband 

satellite"’7.  By  taking  account  of  the  diffuse  auroral  particle 

precipitation  (current)  the  stable  K  x  _R_  diffuse  auroral  geometry 

(corresponding  to  the  observations)  becomes  destabilized  by  this  parallel 

current.  For  a  westward  ambient  d.c.  electric  field,  J^,  and  a  northward 

dominant  electron  density  gradient,  the  relative  drift  velocity  between 

ions  and  electrons  parallel  to  B,  V^,  must  satisfy  the 

condition  -k*v  >k*(cF  /B)  (v. /2. )  for  instability.  The  maximum  growth  rate 
- d  -  -o  i  i,  _,/2 

for  the  instability  is  y  a  n  l(9n  /9y)V,[l  +  SI  v  /Q  v  ]  /2.  The 

o  o  d  e  i  i  e 

instability  is  mainly  field  aligned  (kfl  <<  k^).  The  instability  is  fluid¬ 
like  in  nature  and  so  can  directly  account  for  the  long  wavelength  diffuse 
auroral  scintillation  causing  F  region  irregularities.  In  addition, 
nonlinear  numerical  simulations  of  this  instability,  as  presented  above, 
indicate  equatorward  movement  (convection)  of  plasma  enhancements  and 


Fig.  11  —  Log-log  plots  of  one  dimensional  (a)  x  power  spectrum  P(kx )  and  (b)  y  power 
spectrum  P(ky)  for  L  =  50  km  at  1900  sec.  P(kx )  and  P(ky)  are  obtained  by  averaging 
I  6n  (kx,ky)/N0  I2  over  ky  and  kx  respectively,  where  Sn(x,y)  =  n(x,y)-N0.  In  (a) 
kpx  =  27T/160  km"3  while  in  (b)  kpy  =  2tt/410  km"3 .  The  dots  represent  the  numerical 
simulation  results  while  the  solid  line  is  obtained  from  a  least  squares  fit  (after  ref.  60). 


northward  (poleward)  movement  of  plasma  depletions  (holes).  Furthermore, 
one-dlmenslonal  spatial  power  spectra  of  the  irregularities  in  both  the 

north-south  and  east-west  directions  are  well  described  by  verse  power 

-2 

laws  *  k  ,  for  y  >  4  km  (the  present  wavelength  range  over  which  the 
simulations  have  been  performed). 

Analytic  mode  coupling  theory  has  shown that  the  instability  can 
stabilize  nonlinearlv  by  generating  linearly  damped  harmonics.  With  a 
northward  plasma  density  gradient,  the  dominant  nonlinear  harmonic  is  in 
the  northward  direction;  whereas,  bv  linear  theory  the  dominant  mode  would 

be  in  the  east-west  direction.  The  power  spectral  density  (PSD)  of  the 

_2 

dominant  nonlinear  mode  is  «  k  .  The  nonlinear  evolution  of  two-sided 

(poleward  and  equatorward  plasma  density  gradients)  plasma  enhancements  in 

the  presence  of  field  aligned  currents  and  ambient  auroral  electric  fields 

of  arbitrary  magnitude  and  direction  are  being  studied^.  Results  show 

destabilization  of  the  plasma  slab  on  the  poleward  side  by  a  combination  of 

the  effects  of  convection  (F  x  _B)  and  field  aligned  currents.  Striation- 

like  structures  form  ~  1  hour  with  east-west  PSD  <*  k  n,  n  ■  2-2.5  for  3 

x 

km  <  2ir/k  <  1 00  km  and  north-south  PSD  “  k  ,  m  "  2  for  3 

x  y  ’ 

km  <  2it/k  <  256  km.  In  addition,  the  linear  theory  of  the  CCI  has  been 

So  S9 

extended  *  to  include  other  effects  which  extend  its  range  of 
validity.  Rased  on  all  these  results  the  CCI  appears  to  be  a  relevant 
instability  mechanism  for  high  latitude  scintillation  causing  irregularity 
studies.  Work  is  ongoing  to  extend  the  instability  mechanism  to  regions 
other  than  the  diffuse  auroral  zone.  Numerical  simulations  in  the 
intermediate  wavelength  (100  m-lkm)  regime  are  also  under  way. 

IV.  Conclusions 


In  the  present  paper,  we  have  exhibited  results  from  the  physics 
modelling  of  ionospheric  irregularity  phenomena.  We  have  given  results 
from  three  F  region  ionospheric  areas:  (A)  plasma  cloud  striation 

phenomena,  (B)  equatorial  spread  F  phenomena,  and  (C)  high  latitude  diffuse 
auroral  Irregularity  phenomena.  The  irregularities  in  these  three  areas 
can  produce  deleterious  scintillation  effects  on  satellite  C^T  systems. 
The  three  areas  have  been  studied  by  deriving  the  relevant  plasma  fluid 
equations  with  the  proper  plasma  instability  mechanisms  built  in  to  them 


(physics  modelling)  and  then  studying  these  equations  using  nonlinear 
numerical  simulation  techniques.  In  case  (A)  the  relevant  instability 
mechanism  for  late  time  striation  phenomena  is  the  _E  x  _B  gradient  drift 
instability  mechanism;  in  (B)  it  is  a  form  of  the  collisional  Rayleigh- 
Taylor  instability  mechanism;  and  in  (C)  it  is  the  current  convective 
instability  mechanism  (or  in  the  case  of  the  two-sided  plasma  enhancement, 
a  combination  of  the  E_  x  R_  gradient  drift  instability  and  current 
convective  instability  mechanisms). 

The  studies  presented  show  that  the  theoretical  and  computational 
programs  are  at  a  stage  where  we  can  generate  detailed  information  on  the 
causes  and  structure  of  ionospheric  irregularities.  Meaningful  information 
can  he  obtained  on  realistic  problems.  The  results  provide  a  basis  within 
which  to  interpret  experimental  data  and  act  as  a  vehicle  for  suggesting 
new  experiments.  The  modelling  allows  us  to  attain  a  predictive  capability 
with  respect  to  ionospheric  irregularities  that  can  cause  Cl  satellite 
systems  scinti  llation  phenomena. 
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U.S.  ARMY  BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN  PROVING  GROUND,  MD  2 1005 

OICY  ATTN  TECH  LIBRARY  EDWARD  BAICY 

COMMANDER 

U.S.  ARMY  SATCOM  AGENCY 
FT.  MONMOUTH,  NJ  07703 

OICY  ATTN  DOCUMENT  CONTROL 

COMMANDER 

U.S.  ARMY  MISSILE  INTELLIGENCE  AGENCY 
REDSTONE  ARSENAL,  AL  35809 
OICY  ATTN  JIM  GAMBLE 

DIRECTOR 

U.S.  ARMY  TRADOC  SYSTEMS  ANALYSIS  ACTIVITY 
WHITE  SANDS  MISSILE  RANGE,  NM  88002 
OICY  ATTN  ATAA-SA 
OICY  ATTN  TCC/F .  PAYAN  JR. 

OICY  ATTN  ATTA-TAC  LTC  J.  HESSE 

COMMANDER 

NAVAL  ELECTRONIC  SYSTEMS  COMMAND 
WASHINGTON,  D.C.  20360 

OICY  ATTN  NAVALF.X  034  T.  HUGHES 
OICY  ATTN  PME  117 
OICY  ATTN  PME  I17-T 
OICY  ATTN  CODE  5011 
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COMMANDING  OFFICER 
NAVAL  INTELLIGENCE  SUPPORT  CTR 
4301  SUITLAND  ROAD,  BLDG.  5 
WASHINGTON,  n.C.  20390 

01CY  ATTN  MR.  DUBBIN  STIC  12 

0 ICY  ATTN  NISC-50 

01CY  ATTN  CODE  5404  J.  GALF.T 

COMMANDER 

NAVAL  OCCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  CA  92152 

03CY  ATTN  CODF  532  W.  MOLER 
01CY  ATTN  CODE  0230  C.  BAGGETT 

01CY  ATTN  CODE  81  R.  EASTMAN 

DIRECTOR 

NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 

01CY  ATTN  CODE  4700  TIMOTHY  P.  COFFEY 
26  CYS  IF  UNCLASS.  1  CY  IF  CLASS) 
01CY  ATTN  CODE  4701  JACK  D.  BROWN 
O ICY  ATTN  CODE  4780  BRANCH  HEAD  (150 
CYS  IF  UNCLASS,  1  CY  IF  CLASS) 

0 ICY  ATTN  CODE  7500 

01CY  ATTN  CODE  7550 

0 ICY  ATTN  CODE  7580 

01CY  ATTN  CODE  7551 

0 ICY  ATTN  CODE  7555 

0 [CY  ATTN  CODE  4730  E.  MCLEAN 

0 ICY  ATTN  CODF  4187 

COMMANDER 

NAVAL  SEA  SYSTEMS  COMMAND 
WASHINGTON,  D.C.  20362 
0 ICY  ATTN  CAPT  R.  PITKIN 

COMMANDER 

NAVAL  SPACE  SURVEILLANCE  SYSTEM 
DAHLORFN,  va  22448 

01CY  ATTN  CAPT  J.H.  BURTON 

OFFICER-IN-CHARGE 
NAVAL  SURFACE  WEAPONS  CENTER 
WHITF  OAK,  SILVFR  SPRING,  MD  20910 
0 ICY  ATTN  CODE  F3 1 

DIRECTOR 

STRATEGIC  SYSTEMS  PROJECT  OFFICE 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20376 
0 ICY  ATTN  NSP-2I4I 
niCY  ATTM  NSSP-2722  FRED  WIMBERLY 


COMMANDER 

NAVAL  SURFACE  WEAPONS  CENTER 
DAHLGREN  LABORATORY 
DAHI.GREN,  VA  22448 

01CY  ATTN  COnE  DF-14  R.  BUTLER 

OFFICER  OF  NAVAL  RESEARCH 
ARLINGTON,  VA  22217 
01CY  ATTN  CODE  465 
01CY  ATTN  CODE  461 
01CY  ATTN  CODE  402 
01CY  ATTN  CODE  420 
01CY  ATTN  CODE  421 

COMMANDER 

AEROSPACE  DEFENSE  COMMAND/DC 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  AFB,  CO  80912 

0 ICY  ATTN  DC  MR.  LONG 

COMMANDER 

AEROSPACE  DEFENSE  COMMAND/ XPD 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  AFB,  CO  80912 
0  ICY  ATTN  XPDOO 
01CY  ATTN  XP 

AIR  FORCE  GEOPHYSICS  LABORATORY 
HANSCOM  AFB,  MA  01731 

01CY  ATTN  OPR  HAROLD  GARDNER 
01CY  ATTN  LKB  KENNETH  S.W.  CHAMPION 
01CY  ATTN  OPR  ALVA  T.  STAIR 
0 ICY  A'H'N  PHP  JULES  AARONS 
01CY  ATTN  PHD  JURGEN  BUCHAU 
0  ICY  ATTN  PHD  JOHN  D.  MULLEN 

AF  WEAPONS  LABORATORY 
KIRTLAND  AFT,  NM  87117 
01CY  ATTN  SUL 

0 ICY  ATTN  CA  ARTHUR  H.  GUENTHER 
0 ICY  ATTN  NTYCE  1LT,  G.  KRAJEI 

AFTAC 

PATRICK  AFB,  FL  32925 
0 ICY  ATTN  TF/MAJ  WILEY 
01CY  ATTN  TN 

AIR  FORCE  AVIONICS  LABORATORY 
WRIGHT-PATTERSON  AFB,  OH  45433 
01CY  ATTN  AAD  WADE  HUNT 
0  ICY  ATTN  .AAD  ALLEN  JOHNSON 
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DFPUTY  CHIFF  OF  STAFF 
RESEARCH,  DEVELOPMENT,  &  ACO 
DEPARTMENT  OF  THE  AIR  FORCE 
WASHINGTON,  D.C.  20330 
01CY  ATTN  AFRDO 

HEADOUATERS 

ELECTRONIC  SYSTEMS  DIVISION/XR 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFR,  MA  01731 
01CY  ATTN  XR  J.  DEAS 

HEADOUATERS 

ELECTRONIC  SYSTEMS  DIVISION/YSEA 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01732 
0  ICY  ATTN  YSF.A 

HEADOUATERS 

ELECTRONIC  SYSTEMS  DIVISION/ DC 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFR,  MA  01731 

01CY  ATTN  DCKC  MAJ  J.C.  CLARK 

COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRIGHT- PATTER SON  AFR,  OH  45433 
0 ICY  ATTN  NICD  LIBRARY 

01CV  ATTN  ETDP  R.  BALLARD 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIFFTSS  AFR,  NY  13441 

0 ICY  ATTN  DOC  LIBRARY/TSLD 
0 1 CY  ATTN  OCSE  V.  COYNE 

SAMSO/SZ 

POST  OFFICE  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
(SPACE  DEFENSE  SYSTEMS) 

01CY  ATTN  SZJ 

STRATEGIC  AIR  COMMAND/XPFS 
OFFUTT  AFB,  NB  68113 

0 ICY  ATTN  XPFS  MAJ  B.  STEPHAN 
01CY  ATTN  ADWATF  MAJ  BRUCE  BAUER 
01CY  ATTN  NRT 

01CY  ATTN  DOK  CHIEF  SCIENTIST 

SAMSO/SK 
P.0.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  *>0000 

01CY  ATTN  SKA  (SPACE  COMM  SYSTEMS) 
M.  CLAVIN 


44 


SAMSO/MN 

NORTON  AFB,  CA  92409 
(MTNUTFMAN) 

0 1 CY  ATTN  MNNL  LTC  KENNEDY 

commander 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFR,  MA  01731 

0  ICY  ATTN  EF.P  A.  LORENTZEN 

DEPARTMENT  OF  ENERGY 
LIBRARY  ROOM  G-042 
WASHINGTON,  D.C.  20545 

01CY  ATTN  DOC  CON  FOR  A.  LABOWITZ 

DEPARTMENT  OF  ENERGY 
ALBUOIH'ROUE  OPERATIONS  OFFICE 
P.O.  BOX  5400 
ALBUOUEROUE,  NM  87115 

0 ICY  ATTN  DOC  CON  FOR  D.  SHERWOOD 

EG&C,  INC. 

LOS  ALAMOS  DIVISION 

P.O.  BOX  809 

LOS  ALAMOS,  NM  85544 

0 ICY  ATTN  DOC  CON  FOR  J.  RREEDLOVF. 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 

01CY  ATTN  DOC  CON  FOR  TECH  INFO  DEPT 

01CV  ATTN  DOC  CON  FOR  L-3R9  R.  OTT 

0 ICY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 

01CY  ATTN  DOC  CON  FOR  L-46  F.  SEWARD 

LOS  ALAMOS  SCIENTIFIC  LABORATORY 

P.O.  BOX  1663 

LOS  ALAMOS,  NM  87545 

01CY  ATTN  DOC  CON  FOR  J.  WOLCOTT 

01CY  ATTN  DOC  CON  FOR  R.F.  TASCHEK 

OICV  ATTN  DOC  CON  FOR  F..  JONF.S 

OICY  ATTN  DOC  CON  FOR  J.  MALIK 

01CY  ATTN  DOC  CON  FOR  R.  JEFFRIES 

OICY  ATTN  DOC  CON  FOR  J.  ZINN 

OICY  ATTN  DOC  CON  FOR  P.  KEATON 

OICY  ATTN  DOC  CON  FOR  D.  WESTERVELT 


SANDIA  LABORATORIES 
P.O.  BOX  5800 


ALBUOUEROUE, 

NM 

87115 

OICY 

ATTN 

DOC 

CON 

FOR 

W.  BROWN 

OICY 

ATTN 

DOC 

CON 

FOR 

A.  THORNBROUCH 

OICY 

ATTN 

DOC 

CON 

FOR 

T.  WRIGHT 

OICY 

ATTN 

DOC 

CON 

FOR 

D.  DAHLOREN 

OICY 

ATTN 

DOC 

CON 

FOR 

3141 

OICY 

ATTN 

DOC 

CON 

FOR 

SPACE  PROJECT 

m* 


mm 


mm 


SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

0 ICY  ATTN  DOC  CON  FOR  B.  MITRPHEY 
01CY  ATTN  DOC  CON  FOR  T.  COOK 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  D.C.  20545 

01CY  ATTN  DOC  CON  DR.  YO  SONG 

OTHER  GOVERNMENT 

CENTRAL  INTELLIGENCE  AGENCY 
ATTN  RD/SI,  RM  5G4R,  HO  BLDG 
WASHINGTON,  D.C.  20505 

01CY  ATTN  OSI/PSIP  RM  5F  19 

DEPARTMENT  OF  COMMERCF 
NATIONAL  BITREATT  OF  STANDARDS 
WASHINGTON,  n.c.  20234 

(ALL  CORRES:  ATTN  SEC  OFFICER  FOR) 
01CY  ATTN  R.  MOORE 

INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  &  INFO  ADMIN 
BOULDER,  CO  80303 

01CY  ATTN  A.  JEAN  (UNCLASS  ONLY) 

01CY  ATTN  W.  fJTLAUT 
01CY  ATTN  D.  CROMBIE 
0  ICY  ATTN  L.  BERRY 

NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 
01CY  ATTN  R.  GRUBB 
01CY  ATTN  AERONOMY  LAB  G.  REID 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
OICY  ATTN  I.  GARFUNKEL 
01CY  ATTN  T.  SALMI 
OICY  ATTN  V.  JOSEPHSON 
OICY  ATTN  S.  BOWER 
OICY  ATTN  N.  STOCKWELL 
OICY  ATTN  D.  OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

OICY  ATTN  RADIO  SCIENCES 


BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY,  CA  94701 
OICY  ATTN  J.  WORKMAN 
OICY  ATTN  C.  PRETTIE 

BOEING  COMPANY,  THF 
P.O.  BOX  3707 
SEATTLE,  WA  98124 

OICY  ATTN  G.  KEISTER 
OICY  ATTN  D.  MURRAY 
OICY  ATTN  G.  HALL 
OICY  ATTN  J.  KENNEY 

BROWN  ENGINEERING  COMPANY,  INC. 
CTTMMINGS  RESEARCH  PARK 
HUNTSVILLE,  AL  35807 

OICY  ATTN  ROMFO  A.  DELIBERIS 

CALIFORNIA  AT  SAN  DIEGO,  UNIV  OF 

P.O.  BOX  6049 

SAN  DIEGO,  CA  92106 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  0213H 
OICY  ATTN  D.B.  COX 
OICY  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  20734 
OICY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

OICY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
OICY  ATTN  H.  LOGSTON 
OICY  ATTN  SECURITY  (PAUL  PHILLIPS) 

ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
OICY  ATTN  J.  ROBERTS 
OICY  ATTN  JAMFS  MARSHALL 
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GENERAL  ELECTRIC  COMPANY 
SPACF  DIVISION 
V ALLEY  FORCE  SPACE  CENTER 
CODOARD  BLVD  KINO  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  IS  10 1 

01CY  ATTN  M.H.  BORTNER  SPACE  SCI  LAB 

OF.NERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NY  13201 
0 ICY  ATTN  F.  REIBERT 

OENERAL  ELECTRIC  TECH  SERVICES  CO.,  INC. 
HMES 

COURT  STREET 
SYRACUSF ,  NY  13201 

01CY  ATTN  0.  MILLMAN 

CFNFRAL  RFSEARCH  CORPORATION 
SANTA  BARBARA  DIVISION 
P.O.  BOX  6770 
SANTA  BARRARA,  CA  93111 
0 ICY  ATTN  JOHN  ISF,  JR. 

0 ICY  ATTN  JOEL  OAR BAR I NO 

OEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 

0 ICY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
01CY  ATTN  TECHNICAL  LIBRARY 
0 ICY  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

OTE  SYL VANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

0 ICY  ATTN  MARSHALL  CROSS 

HSS ,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

0  ICY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
ISO  DAVENPORT  HOUSE 
CHAMPAIGN,  TL  61820 

(ALL  COP.RFS  ATTN  DAN  MCCLELLAND) 

0 ICY  ATTN  k.  YEH 


INSTITUTE  FOR  DEFENSE  ANALYSES 
400  ARMY-NAVY  DRIVE 
ARLINGTON,  VA  22202 
0 ICY  ATTN  J.M.  AEIN 
0  ICY  ATTN  ERNEST  BAUER 
01CY  ATTN  HANS  WOLFARD 
01CY  ATTN  JOEL  BENGSTON 

INTL  TEL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

01CY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

11011  TOR RE YANA  ROAD 
P.O.  BOX  85154 
SAN  DIEGO,  CA  92138 

0 ICY  ATTN  J.L.  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
APPLIFD  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LA URAL,  MD  20810 

01CY  ATTN  DOCUMENT  LIBRARIAN 
0 ICY  ATTN  THOMAS  POTEMRA 
01CY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
01CY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O  DRAWER  00) 

SANTA  BARBARA,  CA  93102 
0  ICY  ATTN  DAS I AC 
01CY  ATTN  TIM  STEPHANS 
01CY  ATTN  WARREN  S.  KNAPP 
01CY  ATTN  WILLIAM  MCNAMARA 
01CY  ATTN  R.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
01CY  ATT:  IRWIN  JACOBS 

LOCKHEED  MISSILES  6  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  <>4088 
01CY  ATTN  DEPT  60-12 
0 ICY  ATTN  D.R.  CHURCHILL 
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LOCKHEED  MISSILES  &  SPACE  CO.,  INC. 

3251  HANOVF.R  STREET 
PALO  ALTO,  CA  <14304 

0  ICY  ATTN  MARTIN  WALT  DEPT  52-12 
0 ICY  ATTN  W.L.  IMHOF  DEPT  52-12 
0 ICY  ATTN  RICHARD  0.  JOHNSON  DEPT  52-12 
01CY  ATTN  J.B.  CLADIS  DEPT  52-12 

LOCKHEED  MISSILE  &  SPACE  CO.,  INC. 
HUNTSVILLE  RESEARCH  &  ENGR.  CTR. 

4800  BRADFORD  DRIVE 
HUNTSVILLE,  AL  35807 
ATTN  DALE  H.  DIVIS 


MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

01CY  ATTN  JOHN  MORGANSTF.BN 
01CY  ATTN  G.  HARDING 
0 ICY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
01CY  ATTN  W.  HALL 
01CY  ATTN  W.  FOSTER 


MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  ROX  5837 
ORLANDO,  FL  32805 

0 ICY  ATTN  R.  HEFFNER 


M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

01CY  ATTN  DAVID  M.  TOWLE 
0 ICY  ATTN  P.  WALDRON 
0 ICY  ATTN  L.  LOUOHLIN 
0 ICY  ATTN  0.  CLARK 


MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  02647 
0 ICY  ATTN  N.  HARRIS 
ATTN  J.  *OULE 
ATTN  GEORGE  MROZ 
ATTN  W.  OLSON 
ATTN  R.W,  HALPRIN 
ATTN  TECHNICAL  LIBRARY  SERVICES 


01CY 
0  ICY 
01CY 
01CY 
01CY 


MISSION  RESEARCH  CORPORATION 

735  STATE  STREET 

SANTA  BARBARA,  CA  93101 


PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVERFIELD  BLVD. 

SANTA  MONICA,  CA  90404 
01CY  ATTN  E.C.  FIELD,  JR. 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

01CY  ATTN  IONOSPHERIC  RESEARCH  LAB 


PHOTOMETRICS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02173 

01CY  ATTN  IRVING  L.  KOFSKY 


PHYSICAL  DYNAMICS,  INC. 
P.O.  BOX  3027 
BELLEVUE,  WA  98009 

01CY  ATTN  E.J.  FREMOUW 


PHYSICAL  DYNAMICS,  INC. 
P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 


♦ 

0  ICY 

ATTN 

P.  FISCHER 

R  &  D  ASSOCIATES 

j 

0  ICY 

ATTN 

W.F.  CREVIER 

P.O.  BOX  9695 

i 

01CY 

ATTN 

STEVEN  L.  GUTSCHE 

MARINA 

DEL  REY,  CA  90291 

0  ICY 

ATTN 

D.  SAPPFNFIELD 

01CY 

ATTN 

FORREST  GILMORE 

i 

< 

0  ICY 

ATTN 

R.  BOGUSCH 

01CY 

ATTN 

BRYAN  GABBARD 

«4 

01CY 

ATTN 

R.  HENDRICK 

01CY 

ATTN 

WILLIAM  B.  WRIGHT, 

j 

‘  1 

01CY 

ATTN 

RALPH  KILB 

01CY 

ATTN 

ROBERT  F.  LELEVIER 

nicy 

ATTN 

DAVE  SOWLF. 

01CY 

ATTN 

WILLIAM  J.  KARZAS 

1 

01CY 

ATTN 

F.  FAJEN 

01CY 

ATTN 

H.  ORY 

i 

01  GY 

ATTN 

M.  SCHFIBE 

01CY 

ATTN 

C.  MACDONALD 

1 

i 

,  i 

01CY 

ATTN 

CONRAD  L.  LONGMIRE 

01CY 

ATTN 

R.  TURCO 

01CY 

ATTN 

WARREN  A.  SCHLUETER 

JR. 
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RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
01CY  ATTN  CULLEN  CRAIN 
01CY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

01CY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  FND  AVENUE 
NEW  YORE,  NY  10023 

01CY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 
P.O.  BOX  2351 
LA  JOLLA,  CA  82038 

01CY  ATTN  LEWIS  M.  LINSON 
01CY  ATTN  DANIEL  A.  HAMLIN 
0 ICY  ATTN  E.  FRIEMAN 
0 ICY  ATTN  E.A.  STRAKER 
0  ICY  ATTN  CURTIS  A.  SMITH 
01CY  ATTN  JACK  MCDOUOALL 

SCIENCE  APPLICATIONS,  INC 
1710  OOODRIDGF  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 


0  ICY 

ATTN 

DONALD  MF.ILSON 

01CY 

ATTN 

ALAN  BURNS 

01CY 

ATTN 

G.  SMITH 

0  ICY 

ATTN 

L.L.  COBB 

nicY 

ATTN 

DAVID  A.  JOHNSON 

0  ICY 

ATTN 

WALTER  G.  CHF.SNUT 

0  ICY 

ATTN 

CHARLES  L.  RINO 

nicY 

ATTN 

WALTER  JAYE 

0  ICY 

ATTN 

M.  BARON 

01CY 

ATTN 

RAY  L.  LEADABRAND 

01CY 

ATTN 

G.  CARPENTER 

01CY 

ATTN 

G.  PRICE 

01CY 

ATTN 

J.  PETERSON 

01CY 

ATTN 

R.  HAKE,  JR. 

0  ICY 

ATTN 

V.  GONZALES 

01CY 

ATTN 

D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01CY  ATTN  W.P.  BOOUIST 

TRW  DEFENSE  &  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
01CY  ATTN  R.  K.  PLEBUCH 
01CY  ATTN  S.  ALTSCHULER 
0  ICY  ATTN  D.  DEE 

VISIOYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01R03 
0  ICY  ATTN  W.  REIDY 
01CY  ATTN  J.  CARPENTER 
01CY  ATTN  C.  HUMPHREY 


Code  262 S  -  2J  Copies 


STEWART  RADIANCE  LABORATORY 
UTAH  STATE  UNIVERSITY 
1  nE  ANGELO  DRIVE 
BEDFORD,  MA  O1730 
0 ICY  ATTN  J.  ULWICK 
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IONOSPHERIC  MODELING  DISTRIBUTION  LIST 
(UNCLASSIFIED  ONLY) 


PLFASE  DISTRIRUTF  ONE  COPY  TO  EACH  OF  THE  FOLLOWINO  PEOPLE: 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 
DR.  P.  MANGE  -  CODE  4101 
DR.  R.  MEIER  -  CODE  4141 
DR.  E.  SZUSZCZEWICZ  -  CODE  4187 
DR.  .1.  GOODMAN  -  CODE  4180 
Dr.  R.  RODRIGUEZ  -  CODE  4187 

A.F .  GEOPHYSICS  LABORATORY 
L.G.  HANSCOM  FIELD 
BEDFORD,  MA  0173D 
DR.  T.  ELKINS 
DR.  W.  SWIDFR 
MRS.  R.  SAGALYN 
DR.  J.M.  FORBES 
DR.  T.J.  KENESHEA 
DR.  J.  AARONS 
DR.  H.  CARLSON 
DR.  J.  JASPF.RSF 

CORNELL  UNIVERSITY 
ITHACA,  NY  14850 
DR.  W.F.  SWARTZ 
DR.  R.  SUDAN 
DR.  D.  FARLFY 
DR.  M.  KELLEY 

HARVARD  UNIVERSITY 
HARVARD  SODARE 
CAMBRIDGE,  MA  02138 
DR.  M.r.  McELROY 
DR.  R.  LINDZF.N 

INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY/NAVY  DRIVE 
ARLINGTON,  VA  22202 
DR.  E.  BAUER 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
PLASMA  FUSION  CENTER 
LIBRARY,  NW1 6-262 
CAMBRIDGE,  MA  02139 

NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREFNBELT,  MD  20771 
DR.  S.  CHANDRA 
DR.  K.  MAFDA 
DR.  R.F.  BENSON 


NATIONAL  TECHNICAL  INFORMATION  CENTER 
CAMERON  STATION 
ALEXANDRIA,  VA  22314 
1 2CY  ATTN  TC 

COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 
DR.  T.  CZUBA 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO.  CA  92152 

MR.  R.  ROSE  -  CODE  5321 

NOAA 

DIRECTOR  OF  SPACE  AND  ENVIRONMENTAL 
LABORATORY 
ROULDER,  CO  80302 
DR.  A.  GLENN  JEAN 
DR.  G.W.  ADAMS 
DR.  D.N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.  F.  DONNELLY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  OIIINCY  STREET 
ARLINGTON,  VA  22217 
DR.  H.  MULLANFV 

PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PA  16802 
DR.  J.S.  NISBET 
DR.  P.R.  ROHRBAUGH 
DR.  L.A.  CARPENTER 
DR.  M.  LEE 
DR.  R.  DIVANY 
DR.  P.  BENNETT 
DR.  F.  KLEVANS 

PRINCFTON  UNIVERSITY 
PLASMA  PHYSICS  LABORATORY 
PRINCETON,  NJ  08540 
DR.  F.  PERKINS 
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SCIENCE  APPLICATIONS,  INC.  UTAH  STATE  UNIVERSITY 

1150  PROSPECT  PLAZA  4TH  AND  STH  STREETS 


LA  JOLLA,  CA  92037 

LOGAN, 

UTAH  R4322 

DR.  D.A.  HAMLIN 

DR. 

R.  HARRIS 

DR.  L.  LINSON 

DR. 

K.  BAKER 

DR.  E.  FRIEMAN 

DR. 

R •  SCHUNK 

STANFORD  UNT'/FRSITY 
STANFORD,  CA  94305 
DR.  P.M.  BANKS 


TJ.S.  ARMY  ABERDEEN  RESEARCH 
AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  MD 

DR.  J.  HEIMERL 

UNIVERSITY  OF  CALIFORNIA, 
RERKELEY 

BERKELEY,  CA  94720 
DR.  M.  HUDSON 

UNIVERSITY  OF  CALIFORNIA 
LOS  ALAMOS  SCIENTIFIC  LABORATORY 
J-10,  MS-664 
LOS  ALAMOS,  NM  S7545 
M.  PONGRATZ 
D.  SIMONS 
C.  BARASCH 
L.  DUNCAN 
P.  BERNHARDT 

UNIVERSITY  OF  CALIFORNIA, 

LOS  ANGELES 
405  HILLGARD  AVENUE 
LOS  ANGELFS,  CA  90024 
DP.  F.V.  CORONITI 
DR.  C.  KENNEL 
DR.  A. Y.  WONG 

UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20740 
DR.  K.  PAPADOPOULOS 
DR.  E.  OTT 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 
DR.  N.  ZABUSKY 
DR.  M.  BIONDI 
DR.  E.  OVERMAN 
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